Optical Detection of a Single Nuclear Spin 
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We propose a method to optically detect the spin state of a ^'^P nucleus embedded in a ^®Si matrix. 
The nuclear-electron hyperfine splitting of the '^^P neutral-donor ground state can be resolved via 
a direct frequency discrimination measurement of the ''^P-bound exciton photoluminescence using 
single photon detectors. The measurement time is expected to be shorter than the lifetime of the 
nuclear spin at 4 K and 10 T. 
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The spin qubit embedded in a crystalline host is an at- 
tractive choice for solid-state quantum computation due 
to its long coherence time. The inevitable cost of the spin 
qubit's isolation from the environment is the difficulty of 
measuring its quantum state. Scanning tunnelling mi- 
croscopy and single molecule spectroscopy have pushed 
close to this limit by reporting single nuclear and single 
electron spin resonances , yet measurement of the 

quantum state of a single spin in the solid state remains 
elusive. 

Nuclei in semiconductors, in particular silicon, have 
generated particular interest for solid-state quantum 
computation. The proposal of Kane 0, 0|j in which 
qubits are encoded as the spin state of single embedded 
■^^P nuclei, proposes to solve the single nuclear measure- 
ment problem by adiabatically transferring the nuclear 
state to the spin of the electron bound to the '^^P impu- 
rity and measuring the electron's spin-correlated charge 
state with a single electron transistor (SET). In practice, 
however, charge fluctuations produced by the SET cou- 
ple back to the '^"'^P nucleus via the electron's strong hy- 
perfine interaction, leading to a decoherence source that 
is not present in bulk silicon. Magnetic resonance force 
microscopy has approached single-electron-spin sensitiv- 
ity, but experimental results thus far have shown that 
this probe induces spin relaxation more quickly than the 
needed measurement time 0, casting doubts onpropos- 
als that seek its use for quantum computation [g. The 
use of ensembles of ^®Si nuclei in bulk silicon with no 
metallic gates has been proposed 0, but this scheme 
replaces the measurement problem with the challenge 
of achieving high nuclear polarization in order to main- 
tain scalability. In this Letter, we propose an all-optical 
method to determine the nuclear spin state of an iso- 
lated ^^P impurity in bulk silicon. This method's incor- 
poration into silicon-based quantum computation archi- 
tectures could ease these difficulties. 

In semiconductors, free excitons can be bound to donor 
and acceptor impurity sites, forming bound excitons. 
When a bound exciton (BE) decays radiatively to a neu- 
tral impurity state, its linewidth is characteristically nar- 
row due to the localization of the exciton ^ . If the hyper- 
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FIG. 1: Energy diagram for the neutral donor (P") and its 
bound exciton (P'',X) in a magnetic field. The (P^jX) state 
is populated via capture of a resonantly excited free exciton. 
The ^'^P nuclear state can be determined by the energy dif- 
ference of a and 6. 



fine interaction between the impurity nucleus and either 
the neutral impurity bound carrier or the impurity BE 
is sufficiently strong, the nuclear state of the impurity 
can be determined via the BE photoluminescence (PL) 
energy. The particular case of the '^^P donor in a ^^Si 
matrix is treated below. 

In a magnetic field, the ground state of the '^^P BE, 
notated (P",X), is split into four hole Zeeman levels. The 
ground state of the neutral donor, P*', is split into two 
electron Zeeman levels 0|. These levels are illustrated 
in Fig. n We assume a large applied magnetic field 
such that the clcctron/holc Zeeman interaction is much 
greater than the hyperfine interaction. 

The P'^ state is described by effective mass theory 
(EMT), where the total wavefunction tpc{j^) of the donor 
electron is given by the product of the EMT enve- 
lope function F{r) and the Bloch wavefunction (p{r) — 
M(r)e*'^ "". Both the Bloch and envelope functions are s- 
like, so the dominant part of the hyperfine interaction is 
the Fermi contact term llll. 
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where go is the free-electron g-factor, 7„ is the ^"'^P gy- 
romagnetic ratio, I is the nuclear spin operator, and S 
is the electron spin operator. When an exciton is bound 
to the neutral-donor site, the two electrons form an anti- 
symmetric spin singlet state 0, . Consequently, the 
hyperfine splitting of the BE is determined only by the 
spin of the bound hole. Since the hole Bloch function is 
p-like, the Fermi contact term is negligible, and assum- 
ing that the envelope function is s-likc, the orbital and 
dipolar terms will be much smaller than the P° contact 
hyperfine splitting. Thus, the energy difference between 
the transitions a and b of the (P°,X) PL, shown in Fig. 1, 
is determined entirely by the hyperfine splitting of the P*^ 
state. From Eq. Q, this splitting is 

Ai?HF,po = ^go^^Blnh\A{0)f■ (2) 

The electron density at the ■^^P (r — 0) site has been 
determined via electron spin resonance to be 0.44 x 
10'^'^ cm~^, corresponding to a AE^p po of 60 MHz 
This splitting has also been calculated with reasonable 
accuracy using a corrected envelope function to account 
for the discrepancy in the observed and EMT ionization 
energies [islll^ . 

The (P'',X) state in Si decays primarily via a non- 
radiative Auger process with a lifetime of 300 ns [l^ . 
However, there exists a zero-phonon radiative channel 
with a lifetime of 2 ms • The PL linewidth of a single 
•^^P donor impurity is thus expected to be approximately 
3 MHz, which is much smaller than the hyperfine split- 
ting of 60 MHz. Experimentally, the PL linewidth from 
an ensemble of '^^P impurities was measured to be less 
than 150 MHz at 2 K, which includes an inhomogeneous 
broadening effect and was limited by the spectrometer 
resolution |l9|]. 

After Auger recombination, the electron/hole recombi- 
nation energy is imparted to the second electron, which is 
ionized. To ensure a fast recapture process of the donor 
electron, one can optically excite conduction electrons. 
Simultaneously, one can resonantly excite free excitons 
which will be captured at the neutral-donor site. The 
BE will relax thermally within its lifetime |inl |. and in a 
magnetic field of 10 T at 4 K, the probability for it to oc- 
cupy the lowest Zeeman level is 0.8. Thus, approximately 
400 photons/s are emitted at the desired transitions a 
and b in Fig. ^ 

The extraction of an emitted photon out of the high- 
refractive-index Si substrate as a well-coUimated beam 
for optical detection is not trivial. For this purpose, 
one can incorporate a planar distributed Bragg reflec- 
tor (DBR) cavity, as shown in Fig. Ufa), at the center 
of which a ^^P impurity is embedded. The cavity mod- 
ifies the radiation pattern and concentrates the emitted 
power in the normal direction. With high and low DBR 
refractive indices ni = 3.0 and n2 = 1.5, the output cou- 
pling efficiency (/3-factor) into a beam emitted in a nor- 
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FIG. 2: (a) DBR planar cavity including a '^^P impurity at 
the center of the one-wavelength-thick optical cavity layer, 
(b) 2D photonic crystal structure including a "^^P impurity at 
the central defect 19!|. 

mal direction can be as high as 0.8 for a random dipole 
orientation This efficiency is achieved at the cost 

of a decreased radiative decay rate by a factor of 3 [20j . 
The overall number of PL photons available for optical 
detection would then be 400 x 0.8/3 ~ 100 photons/s. 

Alternatively, a two-dimensional photonic crystal 
structure (Fig. [SJs) can simultaneously enhance the out- 
put coupling efficiency and the radiative decay rate. A 
detailed analysis based on the finite difference time do- 
main (FDTD) method predicts that an optical mode vol- 
ume of 0.5(A/n)'^, a Q value of 30,000, a spontaneous 
emission decay rate enhancement (Purcell) factor of 100, 
and a /^-factor of 1 can be achieved The expected 

photon flux in this case is 4 x 10^ photons/s. 

In order to find the measurement time needed to de- 
termine the nuclear spin state, we calculate the signal-to- 
noise ratio for direct frequency detection. In the scheme 
illustrated by Fig. O the BE PL is collected and sent to 
input a of a Mach-Zender interferometer. The signal is 
split by the first beamsplitter into arms c and d, which 
have a phase difference of lot. Arms c and d are re- 
combined at the second beamsplitter, and the photons in 
outputs e and / are detected by single photon detectors. 
The photocurrents are then subtracted, and the resulting 
current is time-integrated. 

Assuming an input state with mean frequency luq and 
a Lorentzian spectral lineshape with fuU-width-half-max 
7, the average integrated single- photon current is propor- 
tional to 

{I) ^e-''''/'^cos{LUoT). (3) 

The variance of the current is proportional to 

(A/2) = (/2) - {If = 1 - cos2(c.ot). (4) 

The two input frequencies of the BE PL, corresponding 
to the two nuclear spin states, are centered around a 
mean frequency loq; we notate these frequencies uj± = 
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FIG. 3: Mach-Zender interferometer for frequency discrimi- 
nation of BE PL. 



ujQ ± lS.Lo/2, where ftAw — AE'jjp po. The difference in 
the integrated current between the nuclear states is the 
signal, so the power signal-to-noise ratio per photon for 
such a frequency detection scheme is given by 



S _ [sin(iAwT) -sin(-iAwT)]' 
'N ~ cos2(iAwr) -f e-r^ - 1 



(5) 



in which the interferometer is biased such that ujot — 
(m-|-i)7r. The maximum signal-to-noise ratio is 0.084 per 
photon, assuming the upper bound on the PL linewidth 
of 150 MHz and a delay time r = 2 ns. Assuming 100 
photons per second are collected, the minimum integra- 
tion time needed for a signal-to-noise ratio of unity is 
approximately 0.1 s. If we assume an enhanced sponta- 
neous emission decay rate with a 2D photonic structure 
mentioned above and a collection efficiency of 0.5, the 
minimum integration time is reduced to approximately 
10-3 s. 

Any practical detection scheme will suffer from detec- 
tor efficiency and dark count rates. However, negligible 
dark count rates and a detection efficiency of 0.4 can be 
obtained at 1.1 /^m with a superconducting transition 
edge sensor (TES) HH. The effect of finite detector ef- 
ficiency ?7d is to decrease the signal-to- noise ratio by ryd, 
extending the measurement time. 

Sufficient integration of the luminescence signal mea- 
sures the nuclear spin state only if that state is stable 
over a sufficient number of excitation/luminescence cy- 
cles. Previous experiments measuring the resonance of 
single nuclei have shown only distributions over the nu- 
clear states, as in a high-temperature ensemble, due to 
the randomization of the nuclear spins during the inte- 
gration time of the measurement However, in the 
case of 3^P nuclei in silicon, it is known from ESR ex- 
periments that the equilibrium nuclear spin relaxation 
time Ti of isolated '^^P nuclei exceeds 10 hours at 1.25 K 
and 0.8 T ji^. Such a long Ti of solid-state nuclei is 
a combined consequence of the small number of nuclear 
interactions and a small density of states to which the 
nuclear Zeeman energy may be transferred. Cross relax- 
ation with the phosphorous-bound electron, induced by 



the hyperfine interaction with the emission of a phonon, 
is also a slow process. It was measured to be 5 hours at 
1.25 K and 0.8 T and should theoretically scale to 30 sec- 
onds at 4 K and 10 T [iil. Since the P° ground state 
only lasts nanoseconds before each optical re-excitation 
in this scheme, these already long equilibrium relaxation 
times are negligible. The nucleus may only be appre- 
ciably destabilized by the rapid optical excitations we 
introduce. 

We first argue that nuclear spin relaxation due to opti- 
cally excited, delocalized conduction electrons is not ex- 
pected to be a significant effect. The number of conduc- 
tion electrons required for neutralization of the Auger- 
ionized donors in 1 ns is approximately 10^^ cm"'^, as- 
suming a 4 K electron capture cross section of 4 x 
10~^^ cm^ [23. At this temperature, electrons are ex- 
pected to thermalize quickly, in which case the Ti theory 
described in Ref. is applicable. This theory predicts a 
Ti greater than 10^ seconds. Even ^^Si in heavily doped 
(> 10^^ donors/cm^) silicon has a measured Ti in ex- 
cess of 200 minutes at 4 K 0, ; the Ti correspond- 
ing to '^'^P nuclei would only be shorter by a factor of 
729gi/73ip = 4, still leaving this timescale unimportant. 
An argument that Ti due to free excitons should also be 
negligibly long follows a similar reasoning, since only the 
spin of the s-like electron interacts appreciably with the 
nucleus, and Boltzmann statistics may still be assumed. 

Of greater concern is the probability of a nuclear spin 
flip during the capture of a free electron following the 
Auger process. Such a nuclear flip arises in second or- 
der perturbation theory, in which a virtual electron cap- 
ture and a virtual electron/nuclear spin flip-flop in the 
neutral-donor state occur concurrently with the energy 
compensated by the real emitted phonon. To estimate 
the probability of such an event, we note that the en- 
ergy cost of a hyperflne-induced electron/nuclear spin 
flip-flop (approximately the 1.2 meV Zeeman energy of 
the neutral-donor spin) is substantially smaller than the 
donor binding energy (45 meV). We may therefore as- 
sume that the density-of-states factors in Fermi's Golden 
Rule are unchanged between the first and second order 
processes, and that they are independent of the initial 
spin state. We also assume that the optical excitation of 
free carriers is not spin selective, and since Ti for these 
carriers exceeds the capture time, this implies that the 
initial spin polarizations are approximately equal. It fol- 
lows from these assumptions that the probability ratio 
between the first and second order processes is well ap- 
proximated by the ratio between their matrix elements. 
The second order matrix element for an electron/nuclear 
flip-flop process, assuming without loss of generality that 
the nucleus begins in the j-IJ-) state, may be written as 
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in which |P^|-IJ.) describes the ionized donor and |P°t-IJ.) 
describes the neutral donor with electron spin up and 
nucleus spin down. These states are the only impor- 
tant ones for nuclear destabilization, since the |P''iJj-) 
state is unperturbed by Tic- The unperturbed eigenen- 
ergy difference between the intermediate and final state, 
Epo^^ — i<^pO|^, is a sum of the electron Zeeman term, the 
nuclear Zeeman term, and AEnp po; the electron Zeeman 
term dominates this sum. The Hamiltonian term TYn 
refers to the interaction leading to the capture of the free 
electron; thus, a first-order matrix element for a neutral- 
ization process without a flip-flop is written as 



(1) 



P+TJI^P^TJ; 



(7) 



We assume that this Coulombic process has no spin se- 
lectivity. It follows that the probability per transition of 
a flip-flop may be written as 
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This probability can be seen from an alternative view- 
point: if we use second-order time-independent pertur- 
bation theory to calculate the mixing of the |P°T-I|) and 
|P°-Llt) states due to the flip-flop terms of He, and pre- 
sume that capture rates to these perturbed states are the 
same as to the unperturbed states, the same probability 
is obtained. 

The free exciton capture process leads to a similar 
second-order probability for nuclear randomization. In 
this case, however, the binding energy of the free exci- 
ton to the neutral donor is only somewhat larger than 
the electron Zeeman energy, so the density-of-states fac- 
tors in the transition rates can become important. How- 
ever, we do not expect this correction to alter the per- 
transition probability by more than a factor of order 
unity. We thus estimate that at 10 T, the probability 
of a nuclear flip for each free-exciton-capture and free- 
electron-capture process is approximately twice the re- 
sult of Eq. ©I (60 MHz/280 GHz)^ = 5 x 10"^ The 
consequences of this probability will be discussed shortly. 

A nuclear spin flip can also occur in a similar manner 
during a free exciton capture or BE Auger decay due to 
the BE (hole/nuclear) hyperfine interaction. This pro- 
cess is less important because the BE hyperflne interac- 
tion is much weaker. The magnitude of this interaction, 
which includes both dipolar coupling and any contribu- 
tion due to a small s-like component of the hole's Bloch 
wavefunction, may be estimated as 2 MHz from the re- 
sults of muon spin resonance experiments jl7j | . Since the 
BE Zeeman energy is comparable to that of the bound 



electron, the probability of a spin flip due to the BE 
hyperfine coupling is smaller by the approximate factor 
(2 MHz/60 MHz)2 ^ IQ-^. Nuclear spin flip due to the 
P'^ hyperfine coupling during radiative decay has a sim- 
ilar order to that during free electron capture, but since 
radiative decay is 7000 times less frequent, this probabil- 
ity may also be neglected. 

It is unfortunate that the predominant decay mecha- 
nism is non-radiative, since each Auger process increases 
the probability of nuclear randomization without provid- 
ing a signal photon. After 2 x 10^ excitations, at which 
point the probability of nuclear randomization exceeds 
1/10, we can only expect to have collected and detected 
25 photons or 5 x 10"^ photons with a DBR planar cavity 
or a 2D photonic crystal, respectively. This still yields a 
usable signal-to-noise ratio of 2 or 400 for the respective 
geometries. Thus, we expect the measurement-induced 
lifetime of the single nuclear spin to be long enough to 
allow for the measurement of the nuclear spin state. 

In summary, we have proposed utilizing the sharp BE 
PL lines for single nuclear spin measurement. We have 
demonstrated that it is theoretically possible to measure 
a single donor impurity spin in silicon using currently 
available technology - a semiconductor microcavity and 
single-photon detectors. We have analyzed only one type 
of BE in Si, due to the experimentally observed strong 
hyperfine coupling and narrow PL linewidths. However, 
bound excitons exist in abundance in most semiconduc- 
tors. Thus, this technique could be applied to many more 
systems if the crystal quality is sufficiently high and the 
hyperfine coupling and photon flux are sufficiently large. 
We believe this technique holds promise for state readout 
in quantum computers utilizing nuclear qubits in semi- 
conductors. 
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